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Abstract. Resonant spin-flavour (RSF) conversions of supernova neutrinos, which 
is induced by the interaction between the nonzero neutrino magnetic moment and 
supernova magnetic fields, are studied for both normal and inverted mass hierarchy. 
As the case for the pure matter-induced neutrino oscillation (Mikheyev-Smirnov- 
Wolfenstein (MSW) effect), we find that the RSF transitions are strongly dependent on 
the neutrino mass hierarchy as well as the value of 613 . Flavour conversions are solved 
numerically for various neutrino parameter sets, with presupernova profile calculated 
by Woosley and Weaver. In particular, it is very interesting that the RSF-induced 
i/q — !■ D(. transition occurs, if the following conditions are all satisfied: the value of 
fi^B {fii, is the neutrino magnetic moment, and B is the magnetic field strength) 
is sufficiently strong, the neutrino mass hierarchy is inverted, and the value of ^13 
is large enough to induce adiabatic MSW resonance. In this case, the strong peak 
due to original i>c emitted from neutronization burst would exist in time profile of 
the neutrino events detected at the Super-Kamiokande detector. If this peak were 
observed in reality, it would provide fruitful information on the neutrino properties. 
On the other hand, characters of the neutrino spectra are also different between the 
neutrino models, but we find that there remains degeneracy among several models. 
Dependence on presupernova models is also discussed. 



PACS numbers: 95.85.Ry, 13.40.Em, 14.60.Pq, 97.60.Bw 
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1. Introduction 

A core-collapse supernova explosion is one of the most spectacular events in astrophysics, 
and it attracts a great deal of attention from many physicists and astronomers. It 
also produces a number of neutrinos and 99% of its gravitational binding energy is 
transformed to neutrinos. Therefore, neutrinos play an essential role in supernovae, 
and their detection by ground-based large water Cerenkov detectors, such as Super- 
Kamiokande (SK) and Sudbury Neutrino Observatory (SNO), would provide valuable 
information on the nature of neutrinos as well as supernova physics. What we can learn 
from the next galactic supernova has been considered in many articles (for a review, see 
[IJ). For example, we can constrain the properties of neutrino oscillations, such as the 
mixing angle between the first and third mass eigenstates (^13), and the mass hierarchy 
[normal (mi <^ m^) or inverted (mi ^ ma)] |21IS]- 

In addition to the nonzero neutrino masses and mixing angles, the nonzero magnetic 
moment is another nature of neutrinos beyond the standard model of particle physics. 
If neutrinos have a nonzero magnetic moment, it leads to precession between left- 
and right-handed neutrinos in sufficiently strong magnetic fields jH Ej- In general, 
nondiagonal elements of the magnetic moment matrix are possible and neutrinos can be 
changed into different flavours and chiralities [6|. Furthermore, with the additional effect 
of coherent forward scattering by matter, neutrinos can be resonantly converted into 
those with different chiralities [7j by a mechanism similar to the well-known Mikheyev- 
Smirnov-Wolfenstein (MSW) effect [Hj. This resonant spin-flavour (RSF) conversion 
induced by the neutrino magnetic moment in strong magnetic flelds was flrst introduced 
to solve the solar neutrino problem, and actually gave the best flt solution before the 
KamLAND result [Oj. However, the recent KamLAND experiment fO] has shown that 
the large mixing angle (LMA) MSW solution is the most favourable one; the RSF 
mechanism is suppressed at the subdominant level. From the KamLAND negative 
results for the solar antineutrino search, an upper bound on the neutrino magnetic 
moment is obtained, <lx 10~^^/xb, where /xb is the Bohr magneton [TTj. This upper 
bound is comparable to the most stringent limit from the stellar cooling argument, 

/i, <(i-4)xio-ivb m 

Although the RSF mechanism does not work at a dominant level in the Sun, it may 
occur efficiently in a denser environment with stronger magnetic field, which is actually 
expected in the case of core-collapse supernovae. The RSF conversion mechanism in 
supernovae has been investigated by many authors [3 HSl CH CHI. Among them, Ando 
and Sato ^1] have studied the RSF effect using a three-flavour formulation with the 
latest oscillation parameters, and pointed out that the combination of the MSW and 
RSF effects makes the crossing scheme very interesting to investigate. Since the RSF 
conversions are very sensitive to the value of Ye, which is the electron number fraction per 
nucleon, they have also investigated the dependence of the RSF effect on presupernova 
models with solar and zero metallicities pH]- It is concluded that the efficient (either 
complete or incomplete, depending on presupernova models) RSF conversions occur 
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when the supernova magnetic field is sufficiently strong, fi^Bo ^ (10~^^;Ub)(10^° G), 
where Bq is the strength of the magnetic field at the surface of the iron core. 

However, all the past studies of three-fiavour RSF effect were based on the 
assumption that the neutrino mass hierarchy is normal (mi ^ ^3), although the 
inverted mass hierarchy (mi 3> ms) has not been excluded at all. For the pure MSW 
effect, it is well-known that the supernova neutrino signal with the case of inverted 
hierarchy would be very different from that with normal hierarchy [21 IS]- From the 
analogy of the conversion mechanisms between the MSW and RSF effects, it is easily 
inferred that the RSF conversions will be also very sensitive to the mass hierarchy. 
Therefore, in this paper, we study three-ffavour RSF conversions with the inverted 
mass hierarchy using the latest neutrino mixing parameters, and give a comprehensive 
discussion concerning the dependence on the mass hierarchy. Since the RSF conversions 
are extremely sensitive to the presupernova models as noted above, we first adopt the 
model of 15Mq with solar metallicity by Woosley and Weaver jTHI as a reference model; 
the results with the model by Woosley et al jT2j with both solar and zero metallicity 
are addressed later. In particular, we show that the RSF conversion in the case of the 
inverted hierarchy with the large ^13 causes very different neutrino signal from the other 
models, i.e., the appearance of sharp peak of the neutronization burst in the time 
profile. If this case were realized actually, it would be clearly confirmed by not only the 
neutrino spectrum but also the luminosity curve. 

After the completion of our calculation, a paper appeared in which three-fiavour 
RSF effect was studied with both normal and inverted mass hierarchy ^Hj- The authors 
studied rather qualitatively, but they did not obtain the result that the neutronization 
peak in the Ue signal would be detected, when the efficient RSF conversion takes place 
in the case of the inverted mass hierarchy and large di^; our numerical approach as 
well as the qualitative discussion clearly indicate that result. Thus, we stress that the 
point itself is first discovered in the present paper as well as that our study is the first 
numerical one which comprehensively tackles three-fiavour RSF effect with both the 
normal and inverted mass hierarchy. 

This paper is organized as following. In section |2l we give the formulation used 
in our calculation, which includes all three-ffavour neutrinos and antineutrinos. In 
sectionEl models of supernova neutrinos, presupernova structure, and particle properties 
of the neutrino are illustrated. In section HJ we give qualitative discussions concerning 
neutrino conversions in supernova matter both for the normal and inverted mass 
hierarchy, and show results of numerical calculations in sectional Finally, in section IHl 
a simple discussion how to obtain information on the neutrino properties from the 
supernova neutrino signal is presented, and dependence on presupernova models are 
discussed. 
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2. Formulation 

2.1. Interaction with matter and magnetic fields 

The interaction of the magnetic moment of neutrinos and magnetic fields is described 

by 

((^i)R|/^int|(^j)L) = /^ijSx, (1) 

where /iy is the component of the neutrino magnetic moment matrix, B± is the magnetic 
field transverse to the direction of propagation, and (i/)r and are the right- and left- 
handed neutrinos, respectively. If neutrinos arc Dirac particles, right-handed neutrinos 
and left-handed antineutrinos are undetectable (sterile neutrinos), since they do not 
interact with matter. On the other hand, if neutrinos are Majorana particles, z/r are 
identical to antiparticles of z/l and interact with matter. In this paper, we assume that 
neutrinos are Majorana particles. The diagonal magnetic moments are forbidden for 
Majorana neutrinos, and therefore only conversion between different flavours is possible, 

e.g., (i/e)R ^ {i^h,t)l or (l/e)L ^ {i^n,T)R- 

Coherent forward scattering with matter induces an effective potential for neutrinos, 
which is calculated using weak interaction theory. The effective potential due to 
scattering with electrons is given by 

V±± - ±V2Gf (^±^ + 2 sin^ 0^^ n^, (2) 

where Ue is the electron number density, Gp is the Fermi coupling constant, and 
is the Weinberg angle. The ± sign in front refers to u (+) and u (— ) and that in the 
parentheses to Uc (+) and u^,,r {—)■ The difference between e and fi,T neutrinos comes 
from the existence of charged-current interaction. The subscript ±± of V refers to the 
first and the second ± sign. The ordinary MSW effect between Ue and p^^r is caused 
by the potential difference T4 — = y++ — V+- — V^GfUs, while that between 9^ 

and P^,i- hy Ve — Vp,^f — — V = — v^Gpne. To include the RSF effect, which 

causes conversion between neutrinos and antineutrinos, we should take into account the 
neutral-current scattering by nucleons: 

V = V2Gf Q - 2 sin^ - \/2GF^n„ (3) 

where np,n^^ are the proton and neutron number density, respectively. For neutrinos 
we add +V to the potential and for antineutrinos —V. Therefore, the RSF conversion 
between and z/^ ,-, which is important for the case considered in this paper, obeys the 
potential difference 

AV = V,- Vp,^ 

= + v)- - V) 

= - V2G'F-^(l-2ye), (4) 
mN 
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where p is the density, m-^ is the nucleon mass, and 1^ = rie/ (rze + n^) is the number of 
electrons per baryon. (When we obtained equation (0]), we assumed charge neutrahty 
Ue = rip.) 

2.2. The simplest case: ^ v'^ conversion 

In this subsection, we give the simplest formulation between i/g and v'^ with the 
mass squared difference and mixing angle between the first and third mass eigenstate 
Am^3,^i3, and consider the properties of the RSF conversion. Here, v'^ represents mass 
eigenstate in matter which can be obtained by some linear combination of and Vr- 
Although this treatment is obviously insufficient for a realistic discussion, it is still useful 
for an intuitive understanding. First, we qualitatively illustrate how the RSF conversions 
depend on the neutrino mass hierarchy, and then more realistic three-flavour formulation 
is given in the next subsection. The time evolution of the mixed state of and v'^ is 
described by the Schrodinger equation 




— /Uer'-Bj 




(5) 



where r is the radius from the center of the star, /^er' is the transition magnetic moment, 
and AH is defined by 

Ath'^ 

A/7 = cos 2^13 -Ar. (6) 

The resonance occurs when AH = 0; it does not occur if Am\.^ > (normal hierarchy), 
because Ye < 0.5 (therefore, AV^ < 0) is satisfied in stellar envelope. Instead, the 
conversion between z/g and is affected by the resonance since the potential difference 
between them is — V^ ,- = —AV{> 0). On the other hand, in the case of the inverted 
hierarchy, the RSF conversion takes place between and u'^, but not between and z/^. 
This situation is very similar to the case of pure MSW effect; higher MSW resonance 
occurs in antineutrino (neutrino) sector if the mass hierarchy is inverted (normal) [21 El- 
When the resonance is adiabatic (or when the density is slowly changing at that 
point and the magnetic field is strong enough), significant conversion occurs. The 
adiabaticity of the RSF conversion is given by 



7, 



RSF _ (2AicT'5±)^ 



er' 



|dAr/dr| 

where the subscript "res" means that the value is evaluated at the resonance point. 
Therefore, if the magnetic field is sufficiently strong at the resonance point, the z/c ^ z/^ 
conversion occurs completely. 
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2. 3. Three-flavour formulation 

Here, we present the three-flavour (six-component) formulation of neutrino mixing, on 
which our discussions depend: 



V 



dr \ z/ 
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(13) 



IJ- 



(We assume the CP phase 5 = in equation (fT^ for 
simplicity.) 

The resonant flavour conversion basically occurs when the two diagonal elements 
in the matrix in equation (jH} have the same value. There are four relevant resonance 
points, but they depends on the neutrino mass hierarchy as discussed in the previous 



subsection. In the case of the normal mass hierarchy, they are for Uc 



v,, 



(MSW-L) 



(MSW-H), 



v.. 



(RSF-L), and 



(RSF-H). Here, the quantities 



such as z/' 



represent mass eigenstates in matter which can be obtained by linear 



combination of v^iipp) and Vt-{Pt)- This is because the Veipe) state coincides with mass 
eigenstate in matter owing to large matter potential, and the other mass eigenstates 
are obtained by the rotation of Vp,{p^ and Vrivr) basis. Suffixes '-L' and '-H' attached 
to 'MSW and 'RSF' indicate whether the density at the resonance points are lower 
or higher. On the other hand, in the case of the inverted hierarchy, RSF-H occurs for 

whereas the other two resonances occurs for 



v' as well as MSW-H for 



the same conversions. This situation is summarized in table ^ 
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Table 1. Flavour conversions that arc important for each resonance, in both cases of 
the normal and inverted mass hierarchy. 



Resonance 


Normal hierarchy 


Inverted hierarchy 


RSF-H 


Vc ^ v'r 


Vc ^ K 


RSF-L 


Vc ^ v[, 


Vc ^ 


MSW-H 


Vc ^ v'r 


Vc ^ V'^ 


MSW-L 


^ 


^ v[, 



3. Models 

3.1. Original neutrino emission 

We adopt, as original neutrino spectrum as well as luminosity curve, the result of 
numerical simulation by Thompson et al JH]; we use the model calculated for the 
I5M0 progenitor star. Their calculation has particularly focused on shock breakout 
and followed the dynamical evolution of the cores through collapse until the first 250 
ms after bounce. They have incorporated all the relevant neutrino processes such as 
neutrino-nucleon scatterings with nucleon recoil as well as nucleon bremsstrahlung; 
these reactions have recently been recognized to give non-negligible contribution to 
the spectral formation. In figures Q and |2l we show the original luminosity curve and 
number spectrum of neutrinos, respectively. In these figures, z/x represents non-electron 
neutrinos and antineutrinos. 

Neutrino luminosity curve is quite characteristic among different flavours. In 
particular, there is a very sharp peak of called 'neutronization burst', whose duration 
is typically ~ 10 ms and peak luminosity is ~ 10^^ erg/s. This strong peak is illustrated 
as follows. As a supernova shock moves outward, it dissociates nuclei into free nucleons, 
which triggers deleptonization process e~p — >• i/gH; these z/g build up a sea because they 
are trapped and advected with matter. When the shock crosses the neutrinosphere, 
within which the created are trapped, they are abruptly emitted. For the other 
flavours and i^x, there is no such a sudden burst; both luminosities glow rather 
gradually and they are similar to each other. Usually, for both the pure MSW effect and 
the RSF effect with normal hierarchy, the most easily detected flavour are transformed 
from original z/g and [TU [T3] . Thus, luminosity curve does not provide any useful 
information on the flavour conversion mechanism. On the other hand, for the RSF 
effect with inverted mass hierarchy, the conversion z/g ^ — > z/g is considered to occur 
via RSF-H and MSW-H (see table [l}; we expect that this case can be distinguished 
from the luminosity curve. 

The other characteristic that provides information on flavour conversion mechanism 
is hierarchy of the average energy {E^^) < {EpJ < {E^J as clearly seen from figure El 
flavour conversions also cause the spectral exchange. This energy hierarchy is explained 
as follows. Since z/x interact with matter only through the neutral-current interactions 
in supernovae, they are weakly coupled with matter compared to z/g and z/g. Thus the 
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Figure 1. The original luminosity of the emitted neutrinos as a function of time, 
calculated by Thompson et al ^Hl- The progenitor mass is 15Mq. 

neutrinosphere of locates at deeper in the core than that of z/e and z/g, which leads to 
higher temperatures for z/x. The difference between z/g and comes from the fact that 
the core is neutron-rich and z/e couple with matter more strongly, through z/cU e~p 
reaction. 

3.2. Presupernova profiles 

We use the precollapse model of massive stars of Woosley and Weaver [TU1. The 
model is ISMq progenitor star with solar metallicity and it is labeled as W95S. The 
density and profiles are quite important for the flavour conversions because they 
determine the resonance regions as well as whether it is adiabatic or not. We show 
in figure El the |p(l — 21^0) | (responsible for RSF) and pY^ profiles (responsible for 
MSW) of the W95S model. We also show A12 = m^/\ml^cos2ei2/2^/2GYE^ and 
Ai3 = Am^gl cos 26i'i/2\/2G-pE^ as two horizontal bands (the bandwidth comes from 
the energy range 5-70 MeV). At intersections between A12, A13 and |p(l — 2ye)|,pFe, 
the RSF and MSW conversions take place. 

In the previous publication, we have investigated the dependence on presupernova 
models [T3]; in the paper, we have also used recent ISMq presupernova models by 
Woosley et al [T7j with both solar and zero metallicity (W02S and W02Z, where 'S' and 
'Z' denotes solar and zero metallicity, respectively). This is because the RSF conversion 
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Figure 2. Original neutrino spectrum integrated to 0.25 s after core bounce, calculated 
by Thompson et al 19 . The progenitor mass is ISAf©. 



is very sensitive to the deviation of Ye from 0.5 (see equation (jH)), which strongly depends 
on the metalhcities as well as on the weak interaction rates adopted in the simulation 
of stellar evolution. In fact, for the W02 models, in which the authors adopted more 
updated weak rates, the |p(l — 2Fe)| profiles suddenly drop at RSF region, yielding rather 
nonadiabatic conversions; these properties are shown in figure lU Nevertheless, we use 
the W95S model as our reference model because of the reasons listed below. First, since 
the conversion probabilities in RSF-H is highly adiabatic when the magnetic field is 
sufficiently strong (section E)), it is easier for us to investigate the dependence on mass 
hierarchy. Second, although the W02 models were calculated using a recent shell model 
and resulted in substantial revisions to the older data sets in the W95S model, there 
are still many uncertainties concerning nuclear physics; thus, these models cannot be 
considered to be decisive one. Finally, as discussed in the previous paper [12], because 
the lifetime of massive stars which end their life by gravitational collapse is much shorter 
than that of the Sun, the progenitors of the galactic supernovae are considered to be 
younger than the Sun, and therefore, more metal rich. Because the deviation of Y^, from 
0.5 is determined by rarely existent nuclei, the suppression of (1 — 2Yf?) will be weaker 
and the RSF conversion will incline to be more adiabatic, even if W02S model is the 
correct one of the progenitor with solar metallicity. 

Although we use static progenitor models in calculating the flavour transition, in 
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Figure 3. Presupernova profiles (W95S) used in our calculations, which is calculated 
by Woosley and Weaver The density and combination that is responsible 

for the RSF conversions [p{l — Syg)], and that for the MSW conversions (pl^). Two 
horizontal bands represent A12 and A13 (these definitions are given in the text); at 
the intersections between them and the profile curves, the RSF and MSW conversions 
occur. 



fact the density profile changes drastically during a neutrino burst (~ 10 s) owing to 
shock wave propagation. However, our particular interest is within first 0.25 s after 
core bounce because the calculation of original neutrino emission pTOJ ends at that time 
(Kill) . In the previous paper, we have already shown that until 0.25 s, using the static 
presupernova and magnetic field models is considered to be a good approximation [H]; 
this is based on the numerical calculation of the Lawrence Livermore group, which is 
the only group succeeding the shock propagation into the outer envelope. 

3.3. Magnetic fields 

We assume that the global structure of the magnetic field is a dipole moment and the 
field strength is normalized at the surface of the iron core with the values 10^° G (nearly 
complete RSF conversion) as well as G (pure MSW conversion). The reason for this 
normalization is as follows. The magnetic fields should be normalized by fields that 
are static and exist before the core collapse, because those of a nascent neutron star 
can hardly affect the far outer region, where the RSF conversions take place, within 
the short time scale of a neutrino burst. As discussed in the previous subsection, since 
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Figure 4. The same as figure but for (a) W02S and (b) W02Z models. 

the shock wave does not affect the resonance region at < 0.25 s after bounce, it is also 
expected that the magnetic field structure and strength at the resonance points are not 
seriously changed at that time. The strength of such magnetic fields above the surface 
of the iron core may be inferred from observations of the surface of white dwarfs, since 
both are sustained against gravitational collapse by the degenerate pressure of electrons. 
Observations of the magnetic fields in white dwarfs show that the strength spreads in 
a wide range of 10^-10^ G [201 ■ Considering the possibility of the decay of magnetic 
fields in white dwarfs, it is not unnatural to consider magnetic fields up to 10^*^ G at the 
surface of the iron core. Then, in equation (jSJ, 5^ = -Bo(^o/^)^ sin ©, where Bq is the 
strength of the magnetic field at the equator on the iron core surface, tq the radius of 
the iron core, and B the angle between the pole of the magnetic dipole and the direction 
of neutrino propagation. Hereafter, we assume sin0 = 1. 

3.4- Neutrino parameters 

We adopt the realistic neutrino mixing parameters inferred from the recent experimental 
results: for the atmospheric neutrino parameters, jAm^gl = 2.8 x 10~^ eV^,sin^2^23 = 
1.0, and for the solar neutrino parameters, Amfg = 5.0 x 10~^ eV^,tan^6'i2 = 0.42. 

As for still uncertain neutrino properties, we must set some assumptions, e.g., 
whether the mass hierarchy is normal or inverted, as well as is large, enough for the 
MSW-H conversion to be adiabatic, or not. There is also uncertainty concerning the 
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Model 


So [G] 


Mass hierarchy 


sin^ 26*13 


MSW-NOR-S 





Normal 


10-6 


MSW-NOR-L 





Normal 


0.04 


MSW-INV-S 





Inverted 


10-6 


MSW-INV-L 





Inverted 


0.04 


RSF-NOR-S 


IQi" 


Normal 


10-6 


RSF-NOR-L 


IQio 


Normal 


0.04 


RSF-INV-S 


IQi" 


Inverted 


10-6 


RSF-INV-L 




Inverted 


0.04 



neutrino magnetic moment tensor /iy; however, since the only relevant parameter is /iy-B, 
this uncertainty is already included in that of the magnetic field strength (/xy = IQ-^^yLis 
is assumed). Thus, there are 8 parameter sets due to: whether the magnetic field is 
zero (label by 'MSW') or sufficiently strong Bq = 10^° G ('RSF'); the mass hierarchy 
is normal ('NOR') or inverted ('INV'); and sin^ 2^13 = 10"^ ('S') or sin^ 2^13 = 0.04 
('L'). We label one model by connecting these labels using hyphen, e.g., MSW-NOR-S. 
It should be particularly noticed that although we label the model with strong magnetic 
field simply by RSF-, it does not mean that the pure RSF effect occurs; every model 
labeled by RSF- is subject to both the MSW and RSF conversions. We summarize these 
models in table El 

4. Qualitative conversion schemes 

We qualitatively illustrate the conversion scheme for each model. Since Ue is the most 
easily detected fiavour, we focus on the conversion channel that gives Ue appearance at 
the detector. Figure El schematically shows crossings among different mass eigenstates 
in matter, which is helpful for the reader to understand the qualitative discussions in 
this section. 

4.1. Normal mass hierarchy 

First, the case of the normal mass hierarchy is addressed (figure Efa)); this case is 
minutely investigated in [1111121, but we again repeat their discussion. When Bq = 0, 
the produced z/g propagate the supernova envelope without experiencing the MSW 
resonance, and reach the stellar surface as the lightest mass eigenstate z/i. The other 
mass eigenstates U2, t's are originated from z/^ whose fiux is considered to be the same 
at the leading order. Thus, the z/g fiux at the detector is described by 

= |[/el|2F° +(l-|f/einF°, (14) 

where and F are the fiux at production and detection, respectively. Since the 
parameter ^13 is essential for the MSW-H conversion, the results are not sensitive to the 
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Figure 5. Schematic illustration of level crossings for (a) normal and (b) inverted 
mass hierarchy. In this figure, adiabatic conversion means that the neutrinos trace 
the solid curve at each resonance point, while nonadiabatic conversion is shown by the 
dashed line. 



value of 6'i3. 

On the other hand, when Bq = 10^° G, RSF-H becomes almost completely 
adiabatic, whereas RSF-L is highly nonadiabatic. At the RSF-H point, the original 
u'^ are converted into Ue, which propagate as the mass eigenstate owing to the large 
matter potential. These Ug, then, cross the nonadiabatic RSF-L region (which gives no 
essential effects) and escape from the star as Ui. The other mass eigenstates are also 
originated from z/x, thus yielding 

Fp. = Fl (15) 

The value of does not matter also in this case. In consequence, the models MSW- 
NOR-S,L are characterized by equation (fT^ . whereas, RSF-NOR-S,L are characterized 
by equation (fT3j) . 

4.2. Inverted mass hierarchy 

The situation changes dramatically in the case of the inverted mass hierarchy 
(figure Efb)). The pure MSW effect in this case has been studied in detail in the 
literatures [210]; the value of di^ is critical for completeness of the conversions. The z/g, 
which are produced as the lightest mass eigenstate in three antineutrino states, then 
cross the MSW-H resonance region. If this resonance is nonadiabatic (MSW-INV-S), 
the flavour conversion does not take place at the MSW-H point, i.e., Vi. This case 
yields the flux characterized by equation ()14|1 . which has already appeared. On the 
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Table 3. Conversion scheme for each modeL In the second column, the most relevant 
conversion channel, which eventually results in Pi, is given. The flux equation which 
is relevant for each model is summarized in the third column, and then the models are 
categorized into groups labeled by A, B, and C, as shown in the fourth column. 
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other hand for the MSW-INV-L model, MSW-H is completely adiabatic, and therefore 
the conversions Uc i^'r ^ ^'s as well as — > ui^2 occur, which results in 

Fp. = (l-|f/e3p)F°. + |^e3lX. (16) 

Since the value of 6*13 (or It/esP) is strongly constrained to be very small by the reactor 
experiment this expression essentially the same as the previous one 

When Bq = 10^° G (RSF-INV-S, L), the z/e ^ u'^ conversion takes place at the 
RSF-H region. This u'^ further enter the MSW-H region. If this MSW resonance is 
nonadiabatic (-S model), these z/^ becomes at the stellar surface. The other relevant 
conversions are the same as the pure MSW model z/g — ^ z^i, t'^ ^ z^2- Thus, the Ue flux 
is given by 

F,. = \U.i\'Fl + |f/e2pF0 + \U,s\'Fl (17) 

which is almost the same as equation ()14j) . On the other hand, for -L model, the mass 
eigenstates z/^ converted from z/g reach the supernova surface as z^i, because the adiabatic 
MSW-H conversion z/^ ^ z/p takes place. With the remaining channels z/^ — z/2, z/g ^ z/3, 
the observed flux is 

Consequently, the models MSW-INV-S and RSF-INV-S are characterized by equation 
(HH), MSW-INV-L is by (HHD and RSF-INV-L is by (UHl). 

In table El we summarize the relevant conversions stated above for each model. 
Although there are eight models and this rather large number is expected to complicate 
the discussion, each model is consequently categorized into one of three groups A, B or 
C, as indicated in table El Fortunately, this greatly reduces the complexity, but within 
one group, we cannot specify each model by the SK observation. 
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5. Results of numerical calculations 



5.1. Conversion probability 

We calculated equation (jHJ numerically with the models given in section El and 
obtained the conversion probabilities for each flavour. Figure Efa) shows the conversion 
probabilities of original i/^ as a function of radius for the RSF-NOR-L model; the 
neutrino energy is 25 MeV. Because the state is not mass eigenstate and maximally 
mixing with Ur, the probability Pii^fi —>■ v^) should oscillate abruptly. In our calculation, 
however, we have averaged out this local behaviour, since we are only interested in the 
global change of the probabilities which is concerned with resonance. Therefore, the 
survival probability starts from 0.5 rather than 1, because of the maximal mixing with 
v-r- A signiflcant amount of z/^ change into z^e, owing mainly to the RSF-H conversion 
which occurs around ~ O.Oli^Q, and the converted propagates as the mass eigenstate 
due to matter effect, not being disturbed by further resonances, to ~ Q.IRq. Around 
this radius, Ve start to mix with other flavour antineutrinos, reducing the probability 
P{''^ii ~^ ^c)- In flgurelHfb), we show the conversion probabilities of original Pe for the 
RSF-INV-L model. As already discussed in 14. 2| the original z/e are transformed into v^^r 
at RSF-H, and they further change to the favoured flavour z/g iii the MSW-H resonance 
point. 



a) RSF-NOR-L; (b) RSF-INV-L 
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Figure 6. Conversion probabilities as a function of radius for (a) the RSF-NOR-L 
and (b) RSF-INV-L models. The injected neutrino energy is taken to be 25 MeV. 
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We show in figure [7| the same probabihties as those shown in figure IHl but the 
assumed value of 6^13 is small ((a) RSF-NOR-S; (b) RSF-INV-S). Figure [7|^a) indicates 
that the conversion probabilities from the original are almost the same as those for 
RSF-NOR-L model; as we have already noted, the value of ^13 does not matter in the 
case of normal mass hierarchy. On the other hand, the behaviours in figure m^b) is 
substantially different from those in figure IHfb). In particular, the most easily detected 
fiavour is not produced from the original v^- All these characteristics are consistent 
with the simple discussions given in section El 
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Figure 7. The same as figure □ but for (a) the RSF-NOR-S and (b) RSF-INV-S 
models. 



5.2. Neutrino signals at the Super- Kamiokande detector 

With the conversion probabilities given in the previous subsection and the original 
neutrino spectrum by Thompson et al jlH], we calculated the flux of each flavour 
neutrinos on the Earth. (From this point on, we assume the galactic supernova neutrino 
burst and take 10 kpc as a distance to the supernova.) Using this flux and the cross 
section of the relevant neutrino interaction at SK as well as the sensitivity of the detector, 
we can calculate the expected event numbers from a future galactic supernova neutrino 
burst. In this paper, we adopt the most dominant reaction i/gp e+n alone; a cross 
section of the reaction has been calculated in detail [22] • 

Figure |Hfa) shows time evolution of the energy-integrated event as a function of 
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time for each group A, B and C, and we show the same in figure IHl^b) but using equally 
spaced bins. From the time evolution of neutrino events, we cannot discern the groups 
A and B, because they show almost the same time profile. For group C, however, since 
the original are converted into z/g, the time profile shows steep neutronization peak, 
and the event number contained in this peak is expected to be statistically significant 
as clearly seen in figure |Hfb); the event number included in the most prominent three 
bins is ~ 180. If the neutronization peak were actually obtained, it strongly indicates 
that the model group C would be favoured; since the group C contains only one model, 
RSF-INV-L, a great number of problems concerning the neutrino properties would be 
solved at the same time. In that case, the neutrino would have the nonzero magnetic 
moment, the mass hierarchy would be inverted, and the value of ^13 would be large 
enough for MSW-H to be adiabatic. 

The expected event number per unit energy range, which is integrated during the 
first 0.25 s after core bounce, is shown in figures IHl^c) andlH^d), in units of counts/MeV 
and counts/bin, respectively. From these figures, the model group C gives the softest 
spectrum, while B the hardest and A an intermediate one. In addition to the time 
evolution of the neutrino events, the number spectrum would provide useful information 
on the fiavour conversion mechanism. Although the available data are restricted in order 
to avoid uncertainties concerning shock wave propagation, the obtained data would be 
statistically significant. Using the spectrum, degeneracy between the group A and B is 
expected to be broken. 



6. Discussion 



6.1. How far can we probe the neutrino properties from the supernova neutrino 
observation? 

The expected neutrino signal at the SK detector has been investigated thus far. However, 
the mechanism of supernova explosions is quite unclear, since all the reliable numerical 
simulations have not succeeded in pushing the shock wave to penetrate the entire core. 
There may be several unknown processes which we have omitted so far, and the original 
neutrino spectrum as well as its luminosity curve are still controversial. Therefore, we 
cannot trust characteristics of the calculation by Thompson et al ^H] in detail. Instead, 
we use rather simple quantities in order to discuss the conversion mechanisms from the 
neutrino signals; this approach is expected to considerably reduce the dependence on 
supernova models. 

We adopt the following quantities: 

^ Number of events for Ec> 2b MeV 
~ Number of events for < 15 MeV 
Number of events for < t/ms < 75 
Number of events for 75 < t/ms < 150' 
in order to represent the spectral hardness and the peak sharpness of neutronization 
burst, respectively. Since the time of the core collapse would never known with the 
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Figure 8. Neutrino signal at the SK detector, which is evaluated for model groups A, 
B and C: (a) Time evolution of the neutrino signal in unit of [counts/s]; (b) the same 
as (a), but in unit of [counts/bin] with la statistical error bars; (c) number spectrum 
of positrons for first 0.25 s after core bounce in unit of [counts/MeV]; (d) the same as 
(c), but in miit of [comits/bin] with la statistical error bars. 



neutrino signal alone, we take time origin t = when the first neutrino signal is 
detected. The place of the model groups A, B and C on the (-RIk^-^sk) plane are 
shown in figure IHl These groups are well separated from each other, and we expect that 
this particular remains unchanged even if the adopted supernova model is different. 

Although the degeneracy problem within each group cannot be solved by the SK 
observation, which mainly detect Uq, it may be possible if z/g could be detected efficiently. 
SNO is such a detector currently data taking with 1 000 tons of heavy water. The 
supernova z/g can be detected via u^d e~pp reaction. Although the u^. are also 
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Figure 9. The model groups A, B and C, plotted on i?gj^ vs i?gj^ plane. The error 
bars include only statistical errors, and are at Icr level. 



detected through a similar reaction, i/ed — > e+nn, these events could be discriminated 
using delayed coincidence technique. Thus at SNO, we expect that only signal can be 
extracted. Figure HUl is the same as figure El for the SNO detector, but is plotted for the 
models in (a) group A and (b) B. As shown in this figure, the degeneracy problem could 
be solved with the data in principle, but the current smallness of the SNO detector 
prevents a statistically significant discussion. 

6.2. Dependence on presupernova models 

Until this point, we have adopted the presupernova model W95S as our reference model. 
However, as already stated in IH. 21 or investigated in jTHj, the RSF conversions are highly 
dependent on the presupernova profiles. In fact, figures El and |3] show that the relevant 
profile for the RSF effect, p(l — 2Ye), is quite different among these models as well 
as their metallicities. On the other hand, the MSW conversions are insensitive to the 
presupernova models because their relevant profile is pY^ that is not subject to the 
deviation of from 0.5. In particular for the W02 models, the value of p(l — Y^) 
suddenly drops at the RSF-H region, yielding rather nonadiabatic conversions. Since 
RSF becomes incomplete but MSW does not change for these models, the model 
group B, which contains RSF-NOR-S, RSF-NOR-L and MSW-INV-L models, should 
be further divided into two subgroups; we define group B I containing two RSF models, 
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Figure 10. The same as figure El but for the SNO detector. The models are plotted 
withm (a) groups A and (b) B. 



while MSW-INV-L is named group B II. Consequently, we have four characteristic model 
groups A, B I, B II and C; in particular, the groups B II and C contain only one specific 
model each. Within each group, we cannot identify each model from the SK observation. 

Figures ITlT a) and (b) show time evolution of neutrino events at SK, calculated 
with the W02S and W02Z models. Comparing these figures with figure IH^a), the 
neutronization peak is not as sharp as that for the W95S model. This is because 
the RSF-H conversion is not efficient owing to sudden drop of the p(l — 21^) profiles. 
However, since that peak is still prominent and statistically significant, the model group 
C would be discriminated from the other groups using the time profile detected at SK. 

Figures ITlT c) and (d) show the number spectrum of the neutrino signal, which is 
evaluated for the W02S and W02Z models, respectively. The spectra for groups A and 
B II are unchanged from the previous calculation with the W95S model (figure IHl^c)), 
because they are essentially described by the pure MSW effect, which is insensitive to 
presupernova models. Therefore, we focus on the other two groups B I and C. For the 
group B I, the hardness of spectra locates between that of group A and B II, and it 
strongly depends on the adopted metallicities. More detailed discussions concerning 
this group have been already given in our previous paper ^Hl, and we refer the reader 
to the literature. For the group C, the low energy peak disappears from the spectrum 
calculated with the W02S model, refiecting the p{l — 21^) profile of the W02S model 
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Figure 11. Time evolution and number spectrum of neutrino events calculated with 
W02S and W02Z models. 



(figure a)), which significantly reduces the efficiency of the RSF-H conversions. On 
the other hand, this low energy peak remains for the W02Z model, also reflecting the 
p(l — Yc) profile. As clearly shown in figure Efb), it mildly changes in the RSF-H region 
for low energy neutrinos. However, if the energy is beyond some critical value, RSF-H 
becomes completely nonadiabatic, because the p(l — 2Fc) profile abruptly drops to zero. 
Thus, the — > z/g conversion is highly efficient for low energy neutrinos but is highly 
inefficient for high energy ones; for the high energy region the conversion is 

relevant, which leads to double peak profile of the spectrum. 
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7. Conclusion 

In this paper, we investigated the RSF conversions in supernovae for both the normal 
and inverted mass hierarchy. As the case for the pure MSW effect, we found that the 
RSF transitions are strongly dependent on the neutrino mass hierarchy, and also on the 
value of 6*13. We first gave quahtative discussion on the neutrino conversions including 
both the RSF and MSW effect, for eight parameter sets summarized in table |21 From 
that consideration, it was found that these models are categorized into only three groups, 
each of which is expected to show characteristic neutrino signal at the SK detector; we 
named these groups A, B and C, and which parameter set is included in each group is 
summarized in table El 

We, then, presented results of numerical calculations of flavour conversions in 
supernova envelope. The density and profiles of the presupernova star calculated 
by Woosley and Weaver [12] (W95S) was adopted for the calculations. As the magnetic 
field structure, we assumed dipole-type and normalized its strength at the surface of 
the iron core. Using the conversion probabilities calculated by such a procedure and the 
original supernova neutrino spectrum as well as luminosity curve given by Thompson et 
al jini, the expected neutrino signal at the SK detector was estimated. As the result, 
it was found that there are clear difference between the model groups both for the 
spectral shape and the time evolution of the neutrino events. In particular, the model 
group C, which include the RSF-INV-L model alone, shows a sharp neutronization 
peak. Therefore, if this peak were detected in reality from the future galactic supernova 
neutrino burst, it would strongly support the RSF-INV-L model. It would indicate that 
many problems concerning the neutrino property should be solved at the same time, 
i.e., the neutrino have the nonzero magnetic moment, the mass hierarchy is inverted, 
and the value of ^13 is large enough to induce the adiabatic MSW-H resonance. From 
the spectral shape, it would be possible to discern the groups A and B, however, within 
each model group we need other information such as the signal of Vei to discriminate one 
model from the others. Although the SNO detector could detect z/g efficiently, current 
smallness of the detector would prevent the statistically significant discussions. 

We also studied the dependence on presupernova models. As already investigated in 
our previous paper jTH] , the RSF conversion is strongly dependent on the deviation of 1^ 
from 0.5, which is quite sensitive to the metallicities as well as the weak interaction rates 
adopted in the simulation of stellar evolution. The presupernova models by Woosley 
et al jTTj (W02S and W02Z) shows the considerably different p(l — 21^) profile from 
that of the W95S model. For both models, there exists sudden drop of the p(l — 2Yc) 
profile at the RSF-H region, which pushes the RSF conversion rather nonadiabatic, and 
the expected neutrino signal was found to be different from that estimated with the 
W95S model. Still, since the neutronization peak for the RSF-INV-L model exists at 
the statistically significant level, it would keep to be useful method also in these cases. 
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